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Abstract 

The two scale convergence of the solution to a Robin's type-like prob- 
lem of a stationary diffusion problem in a periodically perforated domain 
^ ' is investigated. It is shown that the Robin's problem converges to a prob- 

lem associated to a new operator which is the sum of a standard homoge- 
nized operator plus an extra first order "strange" term; its appearance is 
due to the non-symmetry of the diffusion matrix and to the non rescaled 
^— .• ■ resistivity. 

1 Introduction 

■ Periodic homogenization in perforated media with Robin boundary conditions 
prescribed on the boundary of the holes has been extensively studied by many 
authors and we refer for instance to [7j, [Sj, [T5\, In this paper we study 
the stationary diffusion equation in a periodic perforated body where the heat 

, flow is proportional to the temperature field on the boundary of the holes with 

$-h ' a resistivity having zero average value on the boundary of the reference hole. In 

p5J, the authors studied a model problem for a second-order symmetric elliptic 
operator in a periodically perforated domain with the Robin boundary condition 
prescribed on the boundary of the holes. They use the asymptotic expansion 
technique |T2], HZI to obtain the homogenized problem and they construct 
correctors to justify the expansion and then estimate the error. In this paper, 
we consider a problem but with another configuration, namely the holes may or 
may not be connected and the boundary of the holes may intersect the exterior 
boundary of the body. Moreover we assume that the matrix diffusion of the 
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second-order operator may or may not be symmetric. We use the two-scale 
convergence technique [!]], [EH IH] to obtain the two-scale limit system. 
After the decoupling technique, we show that the homogenized problem contains 
a convective term. Its appearance is due essentially to the general character of 
the matrix diffusion and on the fact that the resistivity function is not rescaled as 
usually assumed when dealing with two-scale convergence on periodic surfaces, 
see for instance [T], g], g], [H)„.. 

The paper is organized as follows: in section [21 we define the geometry of 
the perforated body and we give the Robin boundary-value problem setting. 
Section [3] is aimed at showing the existence and unicity of the solution of the 
Robin problem and obtaining a priori estimates of the solution. The asymptotic 
limit via two-scale convergence procedure is analyzed in section^] We obtain the 
homogenized boundary-value problem which is a second order elliptic operator 
containing first and zero order terms. The latter term is classical see e.g. g]. 
The first order term is a convection one and it is null when the diffusion matrix 
is symmetric or with constant coefficients. 

2 Setting of the Problem 

Let O be a bounded domain in R n of variable x — (xi) 1<i<n (n > 2) with 
a smooth boundary T and e a real parameter taking values in a sequence of 
positive numbers tending to zero. As usual in periodic homogenization let 
Y = [0, 1]™ be the generic unit cell of periodicity in the auxiliary space R™ 
of variable y = (yi) 1<i<n - The cell Y is identified to the unit torus E n /Z n . 
A function defined on R n is said to be y-periodic if it is periodic of period 
1 in each j/j variable with 1 < i < n. In the sequel we will suppose that 
any function defined on Y is extended periodically to the whole space R". If 
E (Z) is a function space (where Z is a subset of Y) we denote E# (Z) := 
{w G E (Z) ; w is extended periodically to R™}. 

Let H, the reference hole, be an open subset of Y with a smooth boundary 
£ and set Y s = Y\cl{H) where cl (•) denotes the closure. Thus Y is partitioned 
as Y = Y s U £ U H . Note that we do not require that H is strictly included in Y. 
As a consequence the periodic extension of H may or may not be connected. Let 
us denote \ (y) the characteristic function of Y s in Y. We define the perforated 
material 

n e = {xen-, x (~) = 1} 

and the one codimensional periodic surface 

E £ = jx S O; j e s| . 

Here Q £ represents the matrix or the solid part of f2, by opposition to the 
holes or the void part that is represented by the open subset H e := f2\cZ (O e ). By 
construction, all of these holes are identical and they are periodically distributed 
in f2 with period e in each Xi-direction. Since we use the two-scale convergence 
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method we do not require that the boundary E e = dH £ does not intersect T. 
As in [2| , We shall use the natural extension by zero of any function defined on 

Let f £ be a given function in L 2 (O e ), g £ be given in L 2 (S e ) such that 

\\fe\\mn e) + V~e\\ge\\L^ E )<C. (1) 

Here and throughout this paper C denotes a positive constant independent of 
e. 

Let A (x. y) — (ay ) 1<i j <n be a real- valued matrix function defined on x Y, 
Y-periodic in the second variable y such that there exists two positive constants 
m, M independent of e satisfying the following inequality: 

m\C\ 2 <(ACX)<M\C\ 2 (2) 

for all C G K n . We suppose that the matrix A lies in C (tyi^ (Y))" . We 
note that no symmetry condition on ^4 is assumed. Let fi(y) € L^ 5 (Y s ) such 
that J y (y) > fiQ > where /xq is independent of e. Let a be a Y-periodic 
measurable bounded function defined on £ such that 

I a(y)da(y)=0. (3) 

Let us decompose the function a into its positive and negative parts as follows: 

a — a + — oT , a + = max (a, 0) , aT — max (—a, 0) . 
Assume that the positive part of a satisfies the condition: 

a + (y) > a > a.e. in E. (4) 
Let us consider the following Robin boundary value problem: 



- div (A e Vu e ) + [i e u £ = f £ in fl £ , (5) 
(A £ \/u £ ) ■ v £ + a £ u £ = eg £ on S e , (6) 
= on T (7) 



where 



A £ (x) = A (x, , fi £ (x) = fx (^j , a £ (x) = a (^j 



and v £ is the unit outward normal to f2 e . 

This problem can be regarded as a simplified model of the condensation of 
stream in a periodic cooling structure (see We can also consider as a model 
for treatment planning hypcrthcmia in microvascular tissue- see, e.g. |1U| . The 
boundary condition JHJ means that the heat flow (A £ "Vu £ ) ■ v £ is proportional 
to the temperature u £ with a periodic resistivity given by the function a £ . In 
many situations, the resistivity function is taken to be e m a £ . Since the operator 
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is of order 2 the interesting cases are then m = —2,-1,0,1 and 2. The case 
to = 2 is trivial since we obtain the classical homogenized equation. This can be 
seen easily by using the asymptotic expansion method. The case m — 1 with 
a e > has been studied by @|, ^21 using the two-scale convergence technique. 
In this situation a e is rescaled since the surface E e is of codimension 1 Here 
we study the case m = 0, i.e. a non rescaled resistivity. We use the same 
technique but with a e changing sign. We show that the assumptions ©, J2J 
and a non-symmetric A £ (x) contribute to the description of the effective thermal 
conductivity with convection. The case m = —1 will be studied in a forthcoming 
paper. Note that the case m = —2 is also trivial since it yields that the effective 
thermal conductivity is 0. 



3 Study of the Problem and A priori Estimates 

Let 

V e = {ve H 1 (O e ) ; v = on T} 
equipped with the scalar product 



Vu (x) Vt> (a;) dx 



and the associated norm ||u||y e = (u,u)y 2 which is equivalent to the .ff^-norm 
thanks to the Poincare inequality. The variational formulation of the boundary- 
value problem ©-0 reads as follows: 



For each e > 0, find u e € V £ such that 
a e (u £ , v) = L £ (v) for any v € V e , 

where a e (•, •) is the bilinear form defined on V e x V e by: 

a e (u,v) = I A £ (x)\7u (x) Vv (x) dx + / \i E (x) u (x) v (x) dx 
a £ (x) u (x) v (x) da e [x] 
and L e (•) is the linear form defined on V e by: 



(8) 



f e (x)v (x) dx + e I g e (x) v (x) dcr e (x) . 



Lemma 1 There exists a positive constant C s independent of e such that for 
every v G V e and for every 6 > we have 



IMIiapj.) < °s \(Se) 1 |M|£ 2( n E ) + (fe) ||Vt>||^ a( n. )r 



(9) 
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Proof. Let us introduce the notation 

v k e (x) =v(e (k + y)) 

where k e K £ = {k £ Z n ; e (Y + k) n ^ 0}. By the change of variable x 
e (k + y) we have 

/ v 2 (x)da £ (x) = [ v 2 (x)da E (x) = e n - 1 J2 f H (v)}* d* (y) 

7 ,/- v o 2-i 



k£K c 



keK e 



From the trace theorem we see that for every S > 



[v k £ (y)} da(y) < C s 



< 



C K 



5- 1 / [v k e (y)} dy + 6 \V y v k (y) \ 2 dy 

JY B + k JY s +k 



5' 1 I v(x) 2 dx + Se 2 I \V x v (x) \ 2 dx 

e(Y s +k) Js{Y s +k) 



Hence 



f v 2 (x)da s (x) < e-^[8-'Yl I 



)e(Y s +k) 



v (x) dx + 



fc 2 V / \V x v(x)\ 2 dx] 



< a 



(Sey L I v 2 (x)dx + (Se) / |Vw (x) \ 2 dx 



The Lemma is proved. 



Lemma 2 Let ^JJi^m > C s ||a||x,oo( S j where C s is the constant given in Lemma 
Then a e (•, •) is coercive on V e . 

Proof. Let v E V £ . Then using we have 

a £ (v,v)>m / | Vu (x) \ 2 dx + /io / v (x) 2 dx — ||a||L°°(s) / v (x) 2 da £ (x) . 
Ja c Jn c Js s 

By JSJ we see that for every 5 > 

Os(U)V) > (m - fcCs||o!||i~(s)) / |Vw (ar) | 2 <iar + 
- (fc) _1 C s ||a|| L oc( S) ^ ^ [v (x)} 2 dx. 

! irri 

Choosing S = — . — . Then 
£ V Mo 



a £ (v,v) > c ||Vu||( i3(nE jj« + |MI|2(n E ) 



where cq is the positive constant given by 

/ C s |ja|j L =o (s) \ 
Co = 1 — mm (to, (jlo) > 0. 

This completes the proof. ■ 

In the sequel we shall assume that the condition ^fji^m > C s ||a||ioo( S ) is 
fulfilled. 

Proposition 3 The variational formulation admits a unique solution u £ £ 
V e . Moreover we have the a priori estimates. 

l|VM e ||(L2(n E ))" + \\Ue\\L 2 (n E ) + l|We||i 2 (S e ) < C. (10) 

Proof. The existence and uniqueness is a straightforward application of Lemma 
^and the Lax-Milgram Lemma. It remains to prove the a priori estimates (llUf) . 
Take v =u e in (JSJ. We have 



/ ( y A E \7u E \7u E + Heii^j dx + a^u 2 e da e (x) 
/ f e u £ dx + / (eg e + a~u e ) u £ da e (x) . 



Let us denote 



A e {u e ) := ||Vu £ ||^ 2(0e)r . + ||u e ||ia ( n e ) + IKIll=(£ e )- 

Then using 10 and (JU we obtain 

A £ (u s )<—( f e u e dx+ (eg e + a~u e ) u e da e (x)) (11) 

where c± = min (to, ao) > 0. Applying Young's inequality on the right hand 
side of ifLTjl. we get 

Ae{u e ) < I [ ^|| /e ||| 2(f2e) + -L|| U£ ||| 2(Qe) 

+ ^\\9e\\h { n s ) + (^2 + IMU-(=)) IKIIi^s.)]. (12) 
But in view of © inequality H12[) becomes 

m*) < ^(^+(^+ini-(-))S)kiiW) 



+ + IMU-(e)) c s ^||v U£ ||^ 2(0e)r ] + a 
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Now, appropriate choice of (3,j,5 yields 

A £ (u e ) = ||Vu E ||^ 2(ne)r + |K||£a (n .) + IKIl!=(E £ ) < C. 

The Proposition is now proved. ■ 

One is led to determine the homogenized problem of Namely we 

study the limiting behavior of the solutions u £ as e tends to zero. This the 
subject of the next section. 

4 Homogenization Procedure 

We shall use the well-known two-scale convergence method that we briefly recall 
here the definition and the main results. 

4.1 Two-scale Convergence 

Definition 4 1. A sequence v £ in L 2 (fl) two-scale converges to vg(x,y) G 
L 2 ((flx Y) and we denote this v £ =4 vq if for any f(x, y) G L 2 (O; C# (Y)), 

lim / v £ {x)ip (x, - ) dx = \ \ v (x,y)(f(x,y)dydx. 
^°Jn \ eJ J n J Y 

2. A sequence v £ in L 2 (S e ) two-scale converges to Vo(x,y) € L 2 ((fl x S) and 
s _ 
we denote this v e =4 «o if for any ip(x,y) G C (fi;C# (Y)j, 

lim / ev £ (x)ip (x, - J da £ (x) = / / vo(x,y)(p(x,y)dcr £ (y) dx. 

Proposition 5 1. For any uniformly bounded sequence v £ in L 2 (fi) one can 
extract a subsequent still denoted by e and a two-scale limit vq G L 2 ((flxY) 
such that v £ =4 Uo- 

2. If v £ is in L 2 (S £ ) such that 

e||«e||i» (E .) < C, 

then one can extract a subsequent still denoted by e and a two-scale limit 

s 

vq G L 2 ((fl x S) such that v £ =4 vq. 

4.2 Two-scale limit system 

By virtue of the estimate QJ and the proposition^ there exists / G L 2 (fl x Y) 
and g G L 2 (ft x E) such that, up to a subsequence, one has 

(x \ S 
-J fe (x) =*x(v)f (x, y) , g £ (x) =3 g (x, y) . (13) 

Furthermore we have . 
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Lemma 6 J^, Let u £ be the solution of 0). Then there exists a 

subsequence still denoted by e and two functions u(x) £ Hq (O), u%(x,y) £ 
L 2 (Cl;H^(Y s )/R) such that 

X (~) Ue (x) =t X (y) u (x) , (14) 

X (|) Vu £ (x) =4 x (tf) (Vu (a:) + V„ux (a;, y)) . (15) 
Moreover we have 

lim / £M e (a;)<y9 (x, — J rfcr e (x) = / / u(x)ip(x,y)da (y) dx (16) 



for every ip £ C(0;C # (F s )). 
Lemma 7 JVe /iaue 



lim 

£^0 



,(x)ip(x)a da s (x) = / / u\(x)ip(x)a (y) da (y) dx 



U £ 



for all ip(x) £ C (fi) . 

Proof. Define a function ip (y) £ Hl(Y s )/M. solution of the problem 

-A0 (y) = 0in F s , 
(V0 (y)) ■ v (y) = a (y) on E, (17) 
2/ 1 — > (y) Y-periodic. 

Such a function exists since satisifies ® which is the compatibility condi- 
tion for the solvability of the problem 1)17(1 . Set if) (y) = Vd and consider the 
function %j) e (x) = if) (-). Then we have 



Vu £ (x)ip £ (x) dx — J u £ (x)t/j £ (x) ■ v (^—^j da £ (x) (18) 

u £ (x)a ^— ^ da £ (x) . 



Passing to the limit in the left hand side of (|18|l and taking into account i|15fl 
we find 



lim / u £ 



,(x)a [ - ) da e (x) = I \ \ (y) (Vu (x) + X7 yUl (x, y)) ip (y) dydx. 
K£/ JnJY 

(19) 

Since u £ Hq (CI) we have 



X (y) Vw (x) ij) (y) dydx = I Vu (x) dx) / \ (y) "0 (y) dy = 
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Hence the right hand side of (|19[) becomes 

X (y) V„ui (x, y) ip (y) dydx. 



IQ, JY 

On the other hand, we have 

X (y) V y ui {x, y) ip (y) dydx 



Q JY 



u\ (x, y) div y 4> (y) dydx 



fi JYi 



Ui (x, y) ip(y)-u (y) da (y) dx 
m (x, y) a (y) da (y) dx 



nJs 



which proves the Lemma. ■ 

Now we are able to give the two-scale limit system: 



Proposition 8 The couple (u,Ui) € Hq (Q) x L 2 \ Q;H^ (Y s ) /RJ is the solu- 
tion of the following two-scale homogenized system : 



-div 



—divy (A (V« + VyUi)) — inflxY s , 
(A(Vu + V v ui) ■ v) +au = onClxH, 
y i — > u\ Y — periodic, 

' x A ( + Vyiii) dyj + Jm + J auida (y) = F 



in f2, 



u = on r 

where p, = J y fx (y) dy and F(x) = J Y X (y) f(%, y) d y + / E g (x, y) da (y) . 



(20) 
(21) 
(22) 

(23) 

(24) 



Proof. Let <p (x) G V (fl) and i Pl (x, y) G V ( fi; (y a )J . Choosing u (jc) 
99 (x) + etpi (a;, as a test function in problem (JSJ, we have 



(x) x ( j) ' (v<p (x) + eV xVl (x, |) + V y (pi (as, J 



+ y (x) u (x) ip (x) + eipi (x, ^ 
-!- / a e (x) u (x) if [x) + e<pi (x, 

x{^) fe{x) ip (x) + etpi (x, '- 



dx 
da e (x) 
dx 



(25) 



g e (x) tp (x) + eyt (x, ~J da £ (x). 
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By virtue of (|15f) the first two terms of the left hand side of l|25l) converges to 
[A (y) [Vu (x) + W yUl (x, y)] [V(p (x) + V y <p! (x, y)] 



+fj, (y) u {x) ip (x)]dydx. (26) 

Taking into account lltj|) . and the lemma[71 the third term of the left hand side 
of tends to 



a(y)u 1 (x,y)tp{x)da(y)dx+ / a (y) u (x) ipi (x, y) da (y) dx. (27) 
Thanks to the right hand side of H25JI converges to 



n Uy, 



f(x,y)dy + / g(x,y)da(y) 



ip (x) dx= F(x)cp (x) dx. (28) 



By the density of X> (ft) x e 2? (ft; (Y s j\ in Hq (ft) x L 2 (ft; (Y s ) /R\ 
we get from the limits (|26f) - H28|) the following two-scale weak formulation system: 



(u, ui) € Hi (ft) x L 2 (ft; (Y s ) /Rj is such that 

In Iy s A I Vu + ^v u i\ I Vu + ^v v i] dydx+ ( 29 ) 
/} J n uvdx + J Q J E auivda (y) dx + J n J E auv\da (y) = J n Fvdx 

for all (v, Vi) G H% (ft) x L 2 (ft; H# (Y s ) /Rj . Integration by parts in J2HJI with 

respect to v\ (v — 0) gives (I20|) - (I22|l and with respect to v (v± — 0) yields 
l|23l) - (|24|l . The Proposition is now proved. ■ 

Thanks to the linearity of the first equation of (12011 we can compute u\(x, y) 
in terms of u (x) as follows: 

n 

u i {x,y) = ^Ck (y) g— (x) + 7(y)u (x) + u (x) (30) 

k=l k 

where for each k the function (j. (y) satisfies the auxiliary problem: 

-div(A(y)V( k (y)) = div(A(y)e k ) in ft x Y s , 
A (y) VCfe (y) ■ v = -A (y) e k ■ v on ft x E, 
2/ — *• Cfe (y) ^-periodic, x e ft. 

where = (<5ifc) 1<i<n , (5i/c is the Kronecker symbol. 
The function j(y) satisfies 

-div (A (y)V 7 (y)) = in X Y s , 
A (y) V7 (y) • v = -a (y) on ft x E, 
?/ ~~ * 7 (y) ^-periodic, x G ft. 
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Finally, inserting the relation l|3()J) into the equation (|23|) yields to the ho- 
mogenized equation 

-div (A hom Vu(j;)) + B ■ Vit(s) + Am (x) = F (x) (31) 

where A hom is the matrix with coefficients 



E 



Y,. 



k=l 

B is the vector with components 



"«(</) i + ^ (y) 



dy 



h = J^a (y) Q (y) da (y) - ^ a ik (y) J^- (y) dy 



a (y) Ci (y) da (y) - \ A (y) ejV7 (y) dy 

A is the real number: 

A = / a (y) 7 (y) da (y) + fi 

A (y) V7 (y) V7 (y) dy + fi 



Thus we have proved the following result 

Theorem 9 Let u e be the solution in V e of the Robin boundary problem f5)) -i^|). 
Then Xs ( x ) u e ( x ) two-scale converges to \ (y) u (%) where u (x) is a solution in 
Hq (Q) of the homogenized problem: 

-div (A hom Vu(x)) + B ■ Vu(x) + Xu (x) = F (x) in O, 
u = on r. 

Remark 10 We observe that the limit equation contains an extra strange 
term of order 1. Namely the convection term B ■ Vit. The vector B depends 
closely on the matrix A and the resistivity function a. For example, if A is 
symmetric then B = 0. Indeed 

[ a (y) Q (y)da = - [ A (y) V 7 (y) VQ (y)dy = - [ A (y) VQ (y) V 7 (y) dy) 
and since A is symmetric 

" (y) Q (y) da = A (y) e t V"/ (y) dy. 
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